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Typeset using REVI^ 



Abstract 



We present the final analysis of the light and strange hadron spectra from 
a full QCD lattice simulation with two degenerate dynamical sea quark 
flavours at /? = 5.6 on a 16^ x 32 lattice. Four sets of sea quark masses 
corresponding to the range .69 < m-^jmp < .83 are investigated. For 
reference we also ran a quenched simulation at P^ff = 6.0, which is the 
point of equal lattice spacing, a~^. 

In the light sector, we find the chiral extrapolation to physical u- and 
d- masses to present a major source of uncertainty, comparable to the 
expected size of unqucnching effects. From linear and quadratic fits wc 
can estimate the errors on the hadron masses made from light quarks to 
be on a 15 % level prior to the continuum extrapolation. For the hadrons 
with strange valence quark content, the Np = 2 approximation to QCD 
appears not to cure the well-known failure of quenched QCD to reproduce 
the physical K — K* splitting. 
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I. INTRODUCTION 



The nonperturbative computation of hadronic properties from Quantum Chromody- 
namics (QCD) presents a major challenge in the unraveling of quark flavourdynamics 
from hadronic experiments. The methods and tools of lattice gauge theory have been 
refined over the past two decades resulting in rather precise results (to the level of a few 
percent accuracy in the physical spectrum of light hadrons, i.e. after chiral and con- 
tinuum extrapolations) within the quenched approximation High statistics quenched 
lattice studies on large lattice volumes revealed that the effects of dynamical fermions 
on spectrum and matrix elements appear to lie within a 10 to 20 % range 

The 'solution' of the full QCD binding problem with lattice methods, on the other 
hand, is still very much lagging behind. This is mainly due to the high cost in compute 
effort to encompass the fermionic determinant in the underlying stochastic sampling 
procedures. The simulation of large lattices in full QCD is definitely a task that requires 
the power of the upcoming teracomputers. Nevertheless, with the computing power 
of some several hundred of teraflops hours it is of considerable interest to tackle QCD 
vacuum polarization effects by looking - on intermediate volumes in the scaling regime 
- at quantities with inherent sea quark dependence such as the vrN a-term, the T^'-mass, 
and the quark spin content of the nucleon. 

A full QCD simulation with Wilson fermions is particularly expensive, as the 
fermionic operator in this case carries more degrees of freedom than in the staggered for- 
mulation, and its chiral extrapolation is more cumbersome as the chiral point fluctuates 
with the gauge field on a finite system. SESAM is a second generation simulation which 
is still exploratory, using Hybrid Monte Carlo (HMC) [Q at /? = 5.6 on 16'^ x 32 lattices. 
For technical reasons, we work with two (degenerate) dynamical fermions, Nf = 2. We 
devised several improvements in order to accelerate the computation of the fermionic 
force [^]. In this way, on the available APEIOO hardware we could achieve HMC 
histories of sufficient lengths for a safe estimate of autocorrelation times. This provides 
a sound basis for the error analysis. 

In a full QCD computation there is no difference between sea quarks, which con- 
tribute to the fermion determinant, and valence quarks, which occur in the hadron 
operators that are employed to excite hadronic states from the QCD vacuum. 

In our Nf = 2 scenario, however, one is forced to introduce 'valence ' quarks different 
from sea quarks^ as soon as one wishes to deal with hadrons carrying strangeness. In 
a recent letter devoted to the determination of the light and strange quark masses, 
we have therefore considered hadronic correlators on a set of three different sea quark 
masses, with valence quark content both equal and different to that of the underlying sea 
quark and presented a consistent approach to analyse such 'semiquenched' data. 

In this paper we will extend that work from three to four different sea quark masses 
and present a detailed study of the light and strange hadron spectra. We shall identify sea 
and valence quarks in the light sector (of u and d quarks) and resort to the semiquenched 
ansatz with respect to the strange quarks, as living in a sea of light quarks. 

For reference, we perform a concomitant quenched simulation on equal lattice spac- 
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ing and volume, at /?e// = 6.00. While unquenching definitely leads to a considerable 
decrease of the light quark mass estimate [0 , we find - within our errors - no visible sea 
quark effects both on the light and strange hadron masses. 



II. SIMULATION DETAILS 
A. Hybrid Monte Carlo 

We have performed a large scale simulation of full QCD at /? = 5.6 with two degener- 
ate flavors of dynamical Wilson fermions. We have generated lattices of extent 16^ x 32 
at four different values of the sea-quark hopping parameter using the ^-version |p of 
the HMC algorithm. The parameters used in the HMC update and the statistics for the 
complete runs on the 256 node APEIOO/Quadrics QH2 are given in table |. 

The CPU costs of the HMC are mostly due to the time consuming repeated solution 
of the linear system M'f M X = $ with M being the Wilson fermion matrix. Throughout 
our simulation we employed the bi-conjugate gradient stabilized algorithm (BiCGStab) 
which has been demonstrated to be the most efficient Krylov sub-space solver for Wilson 
fermion inversions p. Using BiCGStab, we computed the linear system in a two step 
procedure 

M^F = <I>, MX = Y. (1) 

In the first stage of the simulation, we preconditioned by use of the o/e decompo- 
sition of the Wilson fermion matrix M M^. Q, Me = 1 — D^oDoe referred to as 
o/e in table |. In a later stage of the simulation we switched from the thinned o/e 
representation det(Me) to the full fermion determinant det(M) in order to employ the 
locally lexicographic SSOR preconditioner [|10| which has been shown to offer up to a 
factor of 2 less computational costs than o/e preconditioning. In table | we refer to this 
part of the simulation as 'SSOR'. 

As a third improvement of the molecular dynamics part within our HMC, we have 
implemented the chronological start vector guess [|ll|]. The optimal depth of the extrap- 
olation, NcsGi has been determined empirically for each Ksea and with respect to the 
representation of the fermionic determinant as listed in table |. 

We have selected the time step size and the number of molecular dynamics steps, 
Nmdi to yield an acceptance rate of > 70% in the global Monte Carlo decision of HMC. 
With decreasing sea quark mass we can observe a variation of the acceptance rate from 
85% to 73%. We have varied the trajectory length Nmd by numbers uniformly distributed 
in the range ±2^JNmd as recommended in Ref. to avoid deadlocks in periodic orbits 



of phase space due to the presence of well defined Fourier modes. 



This value is at the onset of the (quenched) scaling regime. 
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The chosen stopping accuracy, i?, of the iterative solution of M X = $ is the only 
source of systematic error of the HMC. We have defined the convergence criterion by 
R = ll^^^-"^'! = 10-8 throughout our simulations, working at the level of APE's 32-bit 

precision. Beyond R < 10"^, the difference Hamiltonian AH for the global Monte Carlo 
decision, computed in double precision, does not vary significantly. 

We proceeded adiabatically from large to small sea quark masses and, after ther- 
malizing for more than 500 trajectories at each Ksea, for each sea quark mass, we have 
generated 5000 trajectories. From these correlated samples we have chosen 200 decor- 
related lattices per sea quark mass, see table |. 

B. Error estimates 

Since HMC is a Markov process, one is faced with the problem of autocorrelation of 
the generated series of trajectories. Of course one would like to aim at a decorrelated 
sample of configurations. However, since the generation of trajectories for full QCD is 
extremely costly we cannot afford to skip many trajectories as one can do in quenched 
simulations. In order to control the statistical quality of the measured signals we have 
to carefully study the autocorrelation of the Markov chain. 

We paid attention to keep stable conditions for the HMC dynamics to evolve rather 
than re-tuning HMC parameters during production. This provides the setting for a 
reliable determination of the autocorrelation times related to various hadronic quantities. 

For all four Kgea. values, both exponential and integrated autocorrelation times of 
various gluonic and fermionic observables have been measured. The relevant quantity 
for the error determination is the integrated autocorrelation time Tint- We found Tint, 
which is observable-dependent, to be bound from above by of the smallest eigenvalue 
A of the fermion matrix^, r^^ varies between 15 for Kgea = 0.156 and 30 for Ksea = 0.1575, 
however, the integrated autocorrelation times of most hadronic observables lie well below 
this limit. Therefore, we have decided to analyze every 25th trajectory for spectrum and 
decay constants, after thermalization. 

In order to account for possibly remaining correlations within our hadronic observ- 
ables we have carried out a blocking investigation. For our smallest sea quark mass we 
show in fig.|l| the errors of m^, rrip and rriN, as a function of the blocking size. At block 
size 4 to 6 we find the jackknife errors to run into plateaus. Accordingly, we shall use a 
block size of six throughout our analysis applying the bootstrap procedure. Errors (on 
the blocked data) are obtained from bootstrap samples with 250 entries each. A similar 
analysis of our quenched data shows no increase in error with the block size (quenched 
configurations are generated with an over-relaxed Cabbibo-Marinari heatbath update 
and are separated by 250 sweeps). 



^In ref. we shall present a detailed account of the underlying auto-correlation analysis, 
and we shall propose a scaling rule for the critical slowing down of the HMC. 
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We remark that we have investigated the decorrelation efficiency of the HMC with 
respect to topology on the chosen samples^] [|^. It is gratifying that we could establish 
sufficient tunneling of the topological charge through the topological sectors for the 
four Ksea values investigated. For our smallest quark mass we determined an integrated 
autocorrelation time with respect to the topological charge of Tint ~ 50. Furthermore, 
we analyzed some hadronic quantities — which do not explicitly depend on topological 
effects — according to the topological charge content of the configurations. The result of 
this analysis is that no significant dependence on the topological sector was found. 



C. Hadronic observables 

At each of our four sea quark values we have investigated zero momentum two-point 
functions, 

CAB{t) = j:{o\Ai^)xB{om , (2) 

X 

with hadronic excitation operators x listed in table We combined light-quark 
propagators with hopping parameters equal and different to that of the underlying sea 
quark, thus providing ourselves with fifteen hadronic mass combinations at the two 
heaviest sea quarks and ten at the two lightest (see table |T| for the complete list). 

We use the gauge-invariant Wuppertal-smearing procedure to calculate 

"smeared- local" (st) and "smeared-smeared" (ss) correlators. The smearing parame- 
ter is chosen to be a = 4, with N = 50 iteration steps. In an attempt to further improve 
on our ground state projection we carried out an additional run with 100 smearing iter- 
ations at Ksea = 0.1565; although this rendered a somewhat faster drop into the ground 
state it did not alter our fit results. Plots with Hsca = 0.1565 are from our run with 
= 100. 

Our analyses are based on global masses as extracted from single-exponential ffis to 
the correlators: 

C(tUs = A(e-"^* + e-"^(^-*)) 

C(t)bar = Ae-"^* , (3) 

with T = 32. As a cross check, we also determined effective local masses. For mesons 
they are computed iteratively from the implicit equation 

CAB{t) _ e-™-ffW* + e-™-ffW(^-*) 



CAsit + 1) e-"^c«{t){t+i) + g-meff(t)(r-t-i) 



This investigation is a prerequisite for the investigation of quantities related to topology. 
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while for baryons they are determined in the standard manner from the plateau of local 
masses: 



We use the smeared-smeared data to obtain both masses and amplitudes. The fit ranges 
are determined by keeping the upper limit fixed half-way across the lattice , while the 
lower cut in t is varied in the interval 7 — 10[|. The mass plateau range with the best 
X^/d.o.f- value is selected as fit interval. 

Figure |^ illustrates the quality of our data by showing the different effective local 
masses in comparison to the global masses from correlated fits to the two-point functions, 
in the optimal fit ranges. We find that uncorrelated fits lead to consistent results. 

For future reference our 'raw data' from these mass fits are collected in tables 
to |V111| . By inspection of these tables we retrieve m.„/mp ratios of 0.833(5), 0.809(15), 
0.758(11) and 0.686(11) at K^ea = 0.1560,0.1565,0.1570 and 0.1575 respectively. 

We determine the pseudoscalar and vector decay constants from the respective cur- 
rent matrix elements on the lattice: 



u = —ZAm>), 'f = ^{mi\v). (6) 



The matrix elements are extracted from a direct fit to the ratios 
^..^<»HO|.i,.)e-^^=.^co..(».,^ 

- - 3-V^(0|.i|.)e^-f ^cosKi (,„,,* - f )) , (7) 

where the superscripts I and s denote local and smeared operators, respectively. Note 
that in the second equation the operators VV stand generically for I]fc=i VkVk while the 
operator V on the r.h.s. denotes J2k=i ^k- 

The masses in equations |^ and ^ are fixed to the values obtained from the mass fits 
(given in tables |V| through |V111[ ). 



The 'raw data' for the lattice matrix elements (OIAqItt) and 3~^''^(0|\^'|y) as well as 
f-njZA and l/(/yZy) are collected in tables |I^ to [XI1| . 

The renormalization factors Zj^ and Zy are computed perturbatively, as explained 
in the Appendix. 



III. CHIRAL EXTRAPOLATIONS 



^The smeared-local data yield consistent results but correspond to smaller fit ranges in t. 
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A. The light sector 



We will first present our results for particles and decay constants containing non- 
strange quarks only. It is obvious to identify the degenerate sea quarks in our simulation 
with the u and d quarks since, naively, we expect the lightest sea quarks to make the 
largest effect on the hadronic properties. In this scenario the light hadrons are deter- 
mined from our raw data by a chiral extrapolation in quark mass. We call this setting 
"symmetric" , since it involves data points with equal sea and valence quark masses only. 

At this stage one should remember that full QCD vacuum configurations on differ- 
ent sea quark sectors are manifestly decorrelated. This has some bearing on the error 
analysis of hadron spectra, differently from the quenched situation where one normally 
determines entire hadron mass trajectories configuration wise, with ensuing point to 
point correlations. It will be interesting to trace the impact of this peculiarity on the 
accuracy of hadron masses and decay amplitudes under chiral extrapolation, in the full 
QCD situation. 



1. Masses and decay constants 

The pseudoscalar mass is used to extract the critical hopping parameter k'^^^ while 
the value of the light hopping parameter is set by the condition^: 

"^^'fC'-^ = 0.1785. (8) 

The isospin symmetric bare light quark mass is given by: 

light ^ _ _ /gX 

^ rosea, "'sea 



Linear fits to our data for the pseudoscalar (mpg and vector masses (mv,ss) with 
'^sea = i^vai are shown in fig.H. The resulting parameter values from the extrapolations 



are given in table PCIII| , where we employ the following notations: 

= a + b— (10) 

mv = m""* + cmsea + em^^^ with m^ea = (H) 

We find the pseudoscalar data to be well described by the linear ansatz, the fit yielding 
the critical value of k^p„ to be 



^In this paper we use the convention that physical masses |16| are written in capital letters, 
while lattice masses are denoted by small letters. 
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= 0.158507 (12) 

with x^/d.o.f = 0.6. For the vector particle, both hnear and quadratic parametrizations 
yield acceptable fits, with x^/d.o.f = 0.75 and x^/d.o.f = 0.27, respectively. For the 
light hopping parameter, we find 

K^ght _ 0.158462 (linear) , k^I^' = 0.158471 ttl (quadratic) . (13) 

In the light sector we quote the results from the linear ansatz, using the quadratic fit 
to estimate the systematic uncertainties. This gives the following value for the unrenor- 
malized light quark mass 

flight ^ 0.000901(54) (184) , (14) 

the second error being the systematic uncertainty. Note that this value is consistent 
with our previous estimate, m''^'^* = 0.00088(6) 0], obtained from simulations on three 
sea quark masses. 

We can now predict the nucleon and A masses and the vr and p decay constants by 
chiral extrapolation to the point k,]^^^^. The resulting fit parameters, in the notation of 



of eq. [IT], are collected in table PCII1| . The extrapolations of the baryonic masses are 
visualized in fig.^. It turns out that their rrisea-dependence is by a factor 2 to 3 stronger 
than in the mesonic case, leading to a statistical error on the mass extrapolations for 
nucleon and A of 16 and 22 %, respectively. By comparing the deviations among linear 
and quadratic extrapolations, see table we might estimate a systematic error of 15 
% and 24 %, respectively, which is covered by the statistical error, however. In order to 
put these numbers into perspective, one should be aware that we are extrapolating down 
from m^/mp = .686 on the basis of raw data, which carry statistical errors in the range 
of one to two % (see tables |V| to |V111| ). Just for reference: in state-of-the-art quenched 
simulations [jl| the Tsukuba group achieves statistical errors in the region of .5 to 1 %, 
in the range of rriT^/mp down to a value of 0.4 0]! 

For the decay constants we proceed similarly. The renormalized data are displayed 
in fig. ^ they favour the linear extrapolation (see also table Xlll ). Again, by comparing 



linear and quadratic results, we estimate our systematic uncertainties; they amount to 
15 and 3 % for and fp, respectively. 



B. The strange sector 

So far we have used and ^ to set the values of the hopping parameter values in 
the chiral limit and at the u quark mass. In the following we shall briefiy describe our 
procedure to determine the value of the hopping parameter related to the strange quark 
mass 1^. 

Our simulations are based on two 'active', degenerate sea quarks, which we identify 
with the light quarks. The strange quark in this setting has to be treated as an effectively 
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quenched quark that hves in the sea of the two physical hght quarks. In order to account 
for this situation let us, for the sake of clarity, introduce the generic notation for various 
types of hadron masses appearing in the course of our calculations: 



1. russ - both valence quarks are identical to the sea quark, 

2. rrisv - one valence quark coincides with the sea quark, 

3. rrisjy - both valence quarks differ from the sea quark. 

Note that the 'symmetric extrapolations' operate on the data set rriss and suffice to 
determine both the critical and light hopping parameter values, as discussed in section 

irrrai. 



Let us consider the pseudoscalar masses. In a linear parametrization, the three mass 
types can be written in terms of five slope parameters, c to d": 

^PSjSS ~ ^ n^sea. ; 

"^PS,vv = «'""^sea + &"("^vall + "^val2) • 

In the symmetric situation, nisea = ^vaii = ^vai2, this mapping has to collapse to 
degeneracy on the l.h.s. which leads to constraints on the slopes. As a result one ends 
up with the simple form 

/ m|s 53 \ / a' \ / ^ 
V mis,. J U' - 26' 26' ) V ^vai J 

which can be used for simultaneous fitting in rrisea and m^ai- In the spirit of our approach, 
we will identify the light quark mass with rusea, while the strange quark mass is described 
by Tribal- Note that through the degeneracy requirement, we are effectively left with two 
independent slopes onlyPl, a' and b'. 



Valence quark masses in eq. 15 are defined as: 



"^vai = 7;\- !- I • (16) 



2 '^sca. 



In this setting, with the three types of hadron masses we are in the position to perform 
'semiquenched extrapolations' where valence k, values with /tvai 7^ ^^sea are admitted. 

Within the linear ansatz, other hadronic quantities (like masses and decay constants) 
can be written in terms of mgea and rriyai in the generic form 



^In our previous letter we used three independent such slopes [^. 
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(17) 

We have performed semiquenched fits to eqs.^ and |T7| using the subset of mesonic 

data with Ksea, = ^^vali = ^val2) ^sea = ^^vali 7^ ^^vab ; '^sea 7^ '^vali = l^vah ; described 



by the n combinations in table |lTl[ In order to ensure consistency with the above light 



sector analysis we have used the parameters from the symmetric fits as inputs, namely 

n' — Oh Knight ^Icrit _ ^crit J _ „ 

a — zo, Kgg^ , Kgg^, m — m ,c — c. 

The results of such simultaneous fitting are listed in table pCV| and illustrated in the 
plots of fig. p. All fits are characterized by reasonable x^/d.o.f. As we will discuss 
below, we have tested the stability of the procedure by relaxing the constraints. 

It is obvious how to extend the analysis to nonlinear contributions in rrisea and rriyai'- 

m = m"""* + (c - 2d')mra,,, + 2d'm^^i + {e - f - g')ml^^ + f'm^^^m^^i + g'ml^^ (18) 

The parameters from this nonlinear ansatz, with equal constraints from the above sym- 
metric analysis, and with e' = e, are also included in table PCV| . Note that the coefficients 



of ml^^ and w^sea^^vai both come out to be negative, while the prefactor of m^g^j is positive 
but small. 

Motivated by chiral perturbation theory and quenched QCD one might expect, in- 
stead of eq.|l^, a direct connection on the pseudoscalar mass, according to the form 

m = nf"^ + hml^ , (19) 

which amounts to restricting the parameters 

c' = ha' , d' = bb' . (20) 

To check for the validity of this idea, we have entered our entire data set into a 'scatter- 
plot' with axes my and mpg. Fig. |^ reveals that the entries do not collapse to a single 
line but rather exhibit a clear pattern of sea quark mass dependence, thus ruling out 
the one-slope ansatz eq. |T^. 



1. Determination of K^^'^^'^^^ 

There are three options to fix i^f^^^'^se fj^om the spectrum: 
• from the K* mass by solving 

^V,sv(«;Lf K^t'-'^'^g^) Mk* 



where Kg^a* given by eq. O, 



1.16 , (21) 



^ Obviously this statement can be generalized to any dependence of type m = f{mps) 
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or from the Kaon mass by matching 

3.61 , (22) 



• or from the $ meson mass according to 

"^y2,ss(«sla ) Mp 

It is well known that quenched simulations with Wilson fermions fail to reproduce 
the size of the experimental hyperfine splitting among K and K*. According to the 
results of the CP-PACS collaboration [Q, in the continuum limit and on large lattices, 
the value of Mk* (M$) turns out to deviate by 3% (5%) from experiment when matching 
^strange -j.^ Mj^- On the othcr hand they find Mk* in accord with experiment when using 
M$ as input instead. The deviation is generally attributed to quenching errors. 

In the context of the linear ansatz, the determination of k**''""^'^ from these alternative 



scale choices proceeds directly by explicit use of the fit parameters of table XV. Table 



XVI| lists the resulting values. While the two vector conditions, eqs. pT] and |23, lead to 



consistent results, the K meson mass condition asks for a considerably larger value of 



^strange 



With the numbers for k^^^^^^^s.^ {^om. this table one can proceed to compute the meson 
masses in the strange sector, within the linear ansatz. The results are collected in 



table |XVII| . We find that, contrary to the expectation, the discrepancy between the 
lattice results and the experimental hyperfine splitting remains largely unaltered under 
unquenching. One might be tempted to blame the linear ansatz for this failure. However, 
as can be seen from table PCVI] , the spread in k^^^^^^^s.^ jg means decreased under 

a quadratic extrapolation. We shall come back to this point when we discuss the J- 
parameter. 

If one interprets the spread from the three strange quark mass settings as a systematic 
error, our 'best' value for the strange hopping parameter reads: 

^strange _ Q. 15608(14) (46) , (24) 

which agrees with the mass ratios 

^ = 1.59(17) and ^ = 1.78(22). (25) 
rriK rriK 

The experimental mass ratios are Mk*/Mk = 1.8 and M^/Mk = 2.06. The quoted 
value of ft;*^*'''^"^*' implies a strange quark 

^^ange^2GeV) = 151(30)MeV (26) 
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to be compared to our previous value from the 3 sea quark analysis m^'^^^(2GeV) = 
140(20)MeV. 

The decay constants fx and 1//^ can be determined using the semi-quenched ansatz, 
eq.0, with m replaced by fps and 1/ fv respectively. We compile the results in table 
PCV11| . It turns out that the conditions, eqs. ^ ^ and lead to consistent answers, 



the spread of 3% being well covered by the statistical uncertainty. 



2. Stability of the semiquenched analysis 

By lifting the constraint one can convince oneself in two ways of the stability of the 
light sector, with respect to feedback from the strange sector: 

(i) performing an unconstrained fit to eq. [1^, ond finds (with x^/d.o.f = 22/29) for 
the critical value of K,.^n 



= 0.158497+2, (27) 



which is nicely consistent with the result from the symmetric analysis, eq. [12. 

(ii) An equally satisfying result is achieved with respect to kI,'!]^*: the vector masses 
are reproduced with x^/d.o.f = 17/29, and one obtains the value of the light hopping 
parameter as 



<e* = 0. 158451 (28) 
This number is also in good agreement with the outcome of the symmetric analysis, as 



given in eq. 13 



IV. DISCUSSION 
A. The J-parameter 

In the quenched scenario the dimensionless parameter 

has been proposed as a suitable lattice observable to avoid chiral extrapolations alto- 
gether |1^, on the level of the single mass (i.e. effective k) dependence of the approxi- 



mation. Assuming the vector and pseudoscalar trajectories to be linear one can estimate 
the slope within this expression from experimental mass ratio, such that 

J-^^M^, ^yj^; = 0.48(2) . (30) 
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The quenched lattice value for J generally is of the order of .38, i.e. 25 % below the em- 
pirical estimate. It has generally been surmised that this discrepancy provides evidence 
of a quenching error (see e.g. the review 

Our approach to Nf = 2 full QCD treats light and strange quarks on unequal footing, 
as we associate them with sea and valence quarks, respectively. Basically this induces a 
genuine two-parameter dependence of my on the light and strange quark masses. Thus, 
in the Nf = 2 theory, the notion of one effective hopping parameter is not appropriate 
and the lattice determination of J does not eo ipso enjoy the merit of avoiding chiral 
extrapolation. 



We have demonstrated in section |IIIB 1| that our data does not confirm the single 



slope ansatz, eq.|l^. Thus, in order to avoid the problem of choosing an effective slope 
'^^^ , we calculate the lattice value of J directly from the experimental definition, i.e. 



dM'pg 



we insert our lattice masses, c.f. tabs. pCVll] , prVlll| , on the r.h.s. of eq.|30|. We find 

J = 0.33(3) (31) 
if we define the strange quark mass by eq.^, and 

J = 0.32(3) (32) 



for the condition eq.^. Both values are well below the quenched results. 

We compare this result with the outcome from an analysis restricted to the symmetric 
data. In this case an effective slope value b = ^^^Y can be determined by a linear fit 

to my{rn?pg) on the symmetric data set. The result is h^ym = 1.07(6). Following ref. 
ITtI we set the quark mass by the condition my = 1.8mp5, which corresponds to the 
experimental K* / K mass ratio. This then produces the estimate J=0.40 (2), which is 



significantly above the result of our two-slope analysis, eqs. ^ or |3^. We disfavour this 
approach, however, since the (sea) quark mass, which fulfills the condition my = l.Smps 
on the symmetric line is purely effective and does not correspond to a sea of light u and 
d quarks. 

A third possible way to estimate J in the context of linear extrapolations is to 
apply the above procedure not to the symmetric line, but on each individual line of 
fixed sea quark mass, with subsequent extrapolation of J in Ksea- In this approach, 
however, one has to artificially impose the physical condition my = l.Smps on each 
one of the unphysical sea quark values. One can argue that a possible sea quark effect 
could be easily washed out by such unnaturally guided procedure. And indeed: we 
do not recover any appreciable dependence on the dynamical quark mass with this 
method. Moreover, the numbers, J^Hsea = 0.1560) = 0.34(5), J^Hsea = 0.1565) = 
0.35(4), J{Ksea = 0.1570) = 0.35(5), J{Ksea = 0.1575) = 0.36(5) come out close to the 
quenched values. 

Needless to say, the J-analysis does not provide us with independent information: 
obviously if we had succeeded in predicting the experimental Kaon and K* masses with 
a single value of i^^^^^^^se^ ^i^q results for J, eqs. would agree exactly with J^^^p. 

One might blame the linear ansatz for the failure of J and J^^^f to coincide. In this 
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sense this feature might be considered as an evidence for curvature in the vector particle 
trajectory. In order to explore this possibility, we have carrried out additional quadratic 
fits to the vector particle trajectory both on the symmetric and on the full data set. 
The results for the fit parameters can be found in tables pClll| and pCV| . The coefficients 
of the quadratic terms turn out to be negative, albeit zero within the errors. This then 
even lowers the value of J! 

As yet another alternative we have also used an ansatz with the next to linear order 
in the quark mass ~ m?/'^ . Such a behaviour is expected by chiral perturbation theory 
However, we again find that the J parameter decreases compared to the result of 



the linear ansatz. 



B. Consistency of scale determinations 

The lattice numbers for masses and decay constants can be translated into physical 
results once the lattice cutoff at /5 = 5.6, Np = 2 has been determined. This is done 
by matching the lattice number of one observable with its experimental counterpart. 
Obviously, within a complete numerical solution of QCD, the size of should be 
independent of the choice of the particular observable selected to set the scale. Vice 
versa, a variation of the cutoff with the observable provides another measure for the 
systematic uncertainty of our lattice calculation. 



Table pCVIIl| exhibits the values of the cutoff as obtained by matching physical scales 
inside the light sector: the p mass, the nucleon mass and the pion decay constant. 
Within statistical errors, all results appear consistent. The difference between and 
aj^ however reflects a systematic uncertainty of ~ 20%, this being of course related to 
the the error from the chiral extrapolation, c.f. section pi A| . 

The impact of the uncertainty due to the chiral extrapolation can also be demon- 
strated on itself by making quadratic flts to the vector trajectory (see table pCVlll| ). 



One observes a 10 % change in a~ which goes along with an ampliflcation of the error 
under the quadratic extrapolation from 4 to 15 %. 

The physical predictions are collected in table PCIX| for the light sector, and in table 
X.X| for the masses of particles containing strange quarks. Here we used both from 



linear and quadratic flts in order to test for the stability. It turns out that the admission 
of quadratic contributions in fltting the vector trajectory increases both baryon masses 
by 10 to 15 % and does not reproduce the experimental N — A splitting. The uncertainty 
in the strange sector is clearly dominated by the mismatch of and ^^strange^^*-! 

which is connected with the failure to predict the experimental K — K* splitting. The 
physical results for fx and 1//^ are listed in table PCXI] . As we mentioned above, the 
uncertainty due to the choice of k^^"^^'^^'^ is covered by the statistical errors in this case. 

The analogous quenched results = 6.0, = 2.3 [Gel^], 200 conflgurations of 
16^ X 32 lattices) are contained in table pCXIl| . Notice that the errors on the nucleon and 
A masses are smaller by a factor 2 to 4. 

In flgs. ^ and |^ we present a compilation of the various quantities. We conclude 
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that the data resists to reveal clear sea quark effect on these observables. In particular 
there remains the problem to account for the — A mass splitting. It is unlikely that 
an increase in statistics would remedy the situation. 



V. SUMMARY AND CONCLUSIONS 



We have presented, on moderately sized lattices and at fixed /3, a detailed analysis of 
the light and strange hadron spectra in full QCD. Both meet their particular difficulties: 
while the strange spectrum calculation is hampered by the technical requirement of 
Np = 2, the light baryonic sector faces the problem of considerable variation of hadron 
masses under chiral extrapolation. 

We found that with these limitations we are not able to overcome the well-known 
shortcomings of quenched calculations, namely the underestimation of the K — K* and 

— A spin sphttings. 

The experimental J-parameter cannot be explained in a linear scenario of vector 
trajectories and the admission of higher order terms does not help to improve on the 
situation. In view of this result, it would be highly desirable to make more realistic 
computations by including a third type of active sea quark. 

The issue of the N — A splitting could be considerably clarified by closer approach 
to the chiral limit, on larger lattices. Work along these lines is in progress [|19| . 
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APPENDIX A: RENORMALIZATION CONSTANTS 



We briefly present our method of clioice for the extraction of the renormahzation 
constants ZAy- We use the tadpole improved perturbation theory results from Lepage 
and Mackenzie [Q. The procedure is as follows: 



Use the plaquette values (see table PC1V|) to calculate the value of ay(^^) using: 



1 An 

log(-Spn) = — ay(3.41/a)(l - (1.191 + 0.025nf)av] 
o o 



(Al) 



Convert to AIS scheme using 

1 



1 



a;^(3.41/a) ay(3.41/a) 



+ 0.822. 



(A2) 



Run to a scale -. 

a 

Use tadpole improved perturbation theory: 

Za = 1- 0.31ajjg{-) , Zv 
a 



1 - 0.82a- 



MS\ 



(A3) 



Neglecting the light quark dependence of the plaquette we find 

Za = 0.93 , Zv = 0.82 



(A4) 



for Uf = 2, and Za = 0.94, Zy 
rescale our quark fields: 



0.83 in the quenched case. We also need to 



2k^ ^ - 3fi:/4Kc^. 



(A5) 



Matrix elements in the appendix are listed without the rescaling of the quark fields 
and before applying the renormahzation constants. 
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TABLES 



TABLE I. Simulation parameters and characteristic numbers. 







P = 5.6, Nf = 


2, y X T = 16^ X 32 






0.156 


0.1565 


0.1570 


0.1575 


Algorithm 


o/e 


SSOR 


o/e SSOR 


o/e SSOR 


T 


1 


1 


1 1 


1 0.5 




100 ± 20 100 it 20 


100 it 20 100 ± 20 


100 ±20 71 ± 12 


NcSG 


6 


7 


8 9 


11 3 


# of iter. 


85(3) 


89(6) 


168(5) 125(3) 


317(12) 150(6) 


acc. rate[%] 


85 




84 80 


76 73 


# of traj. 


5000 


5000 


1500 3500 


3000 2000 


# of confs. 


200 


200 


200 


200 


^val - /^val comb. 


15 


15 


10 


10 


TABLE IL The operators studied. 


Mesons 






x\ix)XAiO) 


Pseudoscalar: 






XPsix) = 


= P5 = q'{x)-f^q{x) 


Vector: 






Xvi^) = 


V = q'{x)^''q{x) 


Scalar: 






XSc{x) = 


q'{x)q{x) 


Axial- vector: 






XAx{x) = 


-- Ai" = q'{x)-^57''q{x) 


Baryons 






x\ix)XAiO) 


Nucleon: 






XNix) = 


eabciqaC-f5qb)qc 


A: 






Xa(^) = 


(^abciqaC-f^'qb)qc 


Decay Constants 






xy^),XB(o) 


Pseudoscalar: 






{Al,P,),{PlAo),{Al,Ao) 


Voclor: 
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TABLE III. The run pa,raineters for K^di 





{^val} 


0.156 


{0.156, 0.157, 0.1575, 0.158, 0.1585} 


0.1565 


{0.156, 0.1565, 0.157, 0.1575, 0.158} 


0.157 


{0.1555, 0, 1565, 0.157, 0.1575} 


0.1575 


{0.1555, 0, 1565, 0.157, 0.1575} 



TABLE IV. Integrated autocorrelation times Tint for pseudoscalar, vector and nucleon for 
smeared-focal and smeared-smeared correlators (numbers are in units of HMC time). 



p 




V 


Tint{Mps) 


TintiMv) 


rint{MN) 


B 




sl 


SB 


sl 


ss 


sl 


ss 




5.6 


0.1560 


16^ X 32 


22(9) 


< 25 


< 25 


< 25 


< 25 


< 25 


6 


5.6 


0.1570 


163 X 32 


19(6) 


17(5) 


< 25 


< 25 


< 25 


< 25 


6 


5.6 


0.1575 


16^ X 32 


44(20) 


33(22) 


< 25 


< 25 


37(20) 


32(24) 


7 
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TABLE V. Lattice results for the masses of Pion, Rho, Nucleon and Delta at Kgea = 0.156. 

For all fits wo find 0.1 < /d.o.f < 1. 



nconfigs = 198, nboot = 200, correlated 





(9,15) 




(9,15) 




(9,14) 




(8,14) 










rup 




TUN 




TTIA 




0.1585-0.1585 


0.2937 + 


41 
42 


0.4422 + 


64 
71 


0.676 + 


19 

16 


0.788 + 


21 

19 


0.1585-0.1580 


0.3111 + 


43 
41 


0.4507 + 


62 

67 


0.689 + 


17 

15 


0.796 + 


18 
18 


0.1580-0.1580 


0.3277 + 


41 
40 


0.4591 + 


58 
63 


0.711 + 


16 
13 


0.809 + 


17 
17 


0.1585-0.1575 


0.3279 + 


41 

39 


0.4597 + 


58 
62 


0.702 + 


15 

14 


0.804 + 


18 
17 


0.1580-0.1575 


0.3438 + 


39 
38 


0.4681 + 


60 
56 


0.724 + 


15 
13 


0.817 + 


16 
16 


0.1585-0.1570 


0.3113 + 


39 

38 


0. 1G91 1 


58 
57 


0.715" 


15 
l:^ 


0.812 + 


16 
17 


0.1575-0.1575 


0.3594 + 


36 
36 


0.4771 + 


59 
55 


0.747 + 


14 
12 


0.832 + 


16 
17 


0.1580-0.1570 


0.3595 + 


36 
36 


0.4775 + 


58 
55 


0.737 + 


14 
12 


0.825 + 


16 
16 


0.1575-0.1570 


0.3745 + 


34 
34 


0.4865 + 


59 
54 


0.760 + 


13 
12 


0.840 + 


15 
16 


0.1585-0.1560 


0.3756 + 


36 

35 


0.4884 + 


55 

53 


0.742 + 


13 
13 


0.830 + 


16 
16 


0.1570-0.1570 


0.3892 + 


33 
33 


0.4958 + 


60 

51 


0.783 + 


13 
11 


0.857 + 


14 
16 


0.1580-0.1560 


0.3899 + 


33 
34 


0.4969 + 


57 
51 


0.764 + 


13 
12 


0.843 + 


15 
15 


0.1575-0.1560 


0.4040 + 


33 
31 


0.5059 + 


56 
55 


0.786 + 


12 
12 


0.858 + 


13 
16 


0.1570-0.1560 


0.4179 + 


35 
30 


0.5152 + 


51 
53 


0.808 + 


11 
11 


0.875 + 


14 
15 


0.1560-0.1560 


0.4452 + 


32 
29 


0.5345 + 


54 
49 


0.852 + 


10 
11 


0.910 + 


12 
14 
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TABLE VI. Lattice results for the masses of Pion, Rho, Nucleon and Delta at Kgea = 0.1565. 
Wc find 1.5 < x^/d.o.f < 3 for fits to tt and p , and 0.5 < x^/d.o.f < 1.5 for fits to nucleon and 
A. 



i^scs, = 0.1565 




nconfigs 


= 198, nboot = 200, 


correlated 








(^miiii^max) 


(9,15) 




(9,15) 




(9,14) 




(8,14) 




K1-K2 






rup 




ruN 




rriA 




0.1580-0.1580 


0.3092 1 


61 
50 


0.4408 + 


91 
89 


0.656 + 


19 
21 


0.704 + 


37 
36 


0.1580-0.1575 


0.3257 + 


59 
51 


0.4496 + 


82 

81 


0.673 + 


20 
20 


0.721 + 


35 
34 


0.1575-0.1575 


0.3416 + 


58 
51 


0.4582 + 


81 
81 


0.701 + 


18 
18 


0.748 + 


33 
31 


0.1580-0.1570 


0.3418 + 


58 
52 


0.4588 + 


85 
82 


0.688 + 


19 
20 


0.736 + 


34 
31 


0.1575-0.1570 


0.3572 + 


54 
53 


0.4675 + 


84 
80 


0.716 + 


17 
18 


0.763 + 


31 
26 


0.1580-0.1565 


0.3575 + 


54 
52 


0.4683 + 


86 
84 


0.702 + 


18 

20 


0.750 + 


34 
30 


U.lO (U-U.lO (U 


0.3723 + 


54 
52 


0.4770 + 


81 
77 


n 1 A 1 ~^ 
U. (41 _ 


17 

17 


U. /oy _ 


26 
23 


0.1575-0.1565 


0.3724 + 


54 
53 


0.4772 + 


83 
79 


0.730 + 


17 

18 


0.776 + 


28 
25 


0.1570-0.1565 


0.3871 + 


53 
51 


0.4868 + 


77 
74 


0.755 + 


17 
16 


0.802 + 


26 
23 


0.1580-0.1560 


0.3728 + 


54 
52 


0.4782 + 


87 
82 


0.716 + 


18 
20 


0.763 + 


32 
27 


0.1565-0.1565 


0.4016 + 


52 
50 


0. 19GG t 


70 
(ii) 


0.778;^ 


16 
15 


0.827 + 


23 
22 


0.1575-0.1560 


0.3873 + 


54 
51 


0.4872 + 


79 
76 


0.744 + 


18 
18 


0.789 + 


26 
24 


0.1570-0.1560 


0.4017+ 


52 
50 


0.4968 + 


70 
67 


0.768 + 


17 
17 


0.815 + 


24 
23 


0.1565-0.1560 


0.4157 + 


51 
46 


0.5067 + 


68 
66 


0.792 + 


16 
15 


0.839 + 


23 
21 


0.1560-0.1560 


0.4295 + 


50 
42 


0.5168 + 


63 
63 


0.814 + 


15 
15 


0.862 + 


20 
20 
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TABLE VII. Lattice results for the masses of Pion, Rho, Nucleon and Delta at Kgea = 0.157 
We find 1.5 < x^/d-o.f < 3 for fits to p and A , and 0.5 < x^/d.o.f < 1.5 for fits to tt and 



nucleon. 



''^sca — 0.157 




nconfigs 


= 198, nboot = 200, 


correlated 








(^miin^max) 


(8,15) 




(8,15) 




(9,14) 




(9,14) 




K1-K2 






nip 




ruN 




rriA 




0.1575-0.1575 


0.3163 + 


53 
45 


0.4397 + 


72 
76 


0.695 + 


15 
19 


0.774 + 


13 
14 


n 1 K'TK n 1 K'7r\ 
U.iD/o-U.io/U 


0.3328^ 


50 
42 


0.4499 + 


73 
74 


U. / Uo 


15 
15 


U. / OD _ 


12 
14 


0.1570-0.1570 


0.3486^ 


49 
41 


0.4600 + 


70 
72 


0.732 + 


14 
14 


0.805 + 


13 
13 


0.1575-0.1565 


0.3489 + 


51 
40 


0.4603 + 


72 
76 


0.722 + 


15 
15 


0.798 + 


12 
14 


0.1570-0.1565 


0.3641 + 


47 
37 


0.4703 + 


67 
69 


0.746 + 


14 
13 


0.816 + 


14 
13 


0.1565-0.1565 


0.3790 + 


46 
38 


0.4805 + 


64 
68 


0.769 + 


12 

12 


0.835 + 


14 
13 


0.1575-0.1555 


0.3797 + 


48 
37 


0.4812 + 


66 
72 


0.750 + 


14 

13 


0.819 + 


14 
14 


0.1570-0.1555 


0.3939 + 


46 
36 


0.4911 + 


61 
72 


0.774 + 


13 
12 


0.837 + 


14 
13 


0.1565-0.1555 


0.4080 + 


45 
34 


0.5010 + 


63 
73 


0.797 + 


12 
11 


0.856 + 


13 
13 


0.1555-0.1555 


0.4354 + 


42 
30 


0.5211 + 


61 
70 


0.842 + 


10 
10 


0.895 + 


14 
13 
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TABLE VIII. Lattice results for the masses of Pion, Rho, Nucleon and Delta at 



K«,a = 0.1575. For all fits wo find 0.5 < x'^/(\.Q.i< 1.5. 



Ksea = 0.1575 




nconfigs 


= 198, nboot = 


200, correlated 








(+ + \ 
i^''min)''maxj 


(7,15) 




(10,15) 


(8,14) 




(8,14) 














rriN 




771 A 




0.1575-0.1575 


0.2803 + 


45 
23 


U.4Uo ( _ 


55 
60 


0.633 + 


12 

12 


0.695 + 


16 

16 


0.1575-0.1570 


0.2986 + 


42 
20 


0.4178 + 


56 

50 


0.648 + 


10 

11 


0.707 + 


15 
16 


0.1570-0.1570 


0.3159 + 


37 
20 


0.4272 + 


52 
50 


0.673 + 


8 
9 


0.730 + 


14 
14 


0.1575-0.1565 


0.3160 + 


37 
19 


0.4275 + 


54 
51 


0.662 + 


9 
10 


0.719 + 


15 
15 


0.1570-0.1565 


0.3325 + 


33 
18 


0.4372 + 


50 
51 


0.687 + 


8 
9 


0.742 + 


13 
14 


0.1565-0.1565 


0.3185 + 


30 
20 


0.1 172 " 


52 
17 


0.712;^ 


6 


0.7G1 + 


12 
13 


0.1575-0.1555 


0.3489 + 


31 
19 


0.4479 + 


55 
48 


0.690 + 


8 
9 


0.743 + 


14 
15 


0.1570-0.1555 


0.3642 + 


30 
20 


0.4577 + 


52 
48 


0.715 + 


7 

8 


0.765 + 


13 
13 


0.1565-0.1555 


0.3790 + 


28 
22 


0.4677 + 


56 
47 


0.738 + 


7 
7 


0.786 + 


11 
13 


0.1555-0.1555 


0.4079 + 


27 
23 


0.4881 + 


47 
43 


0.785 + 


6 
7 


0.828 + 


10 
12 
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TABLE IX. Lattice results pseudoscalar and vector meson decay constants at Kgea = 0.156. 

For all fits wo find 1 < x^/d.ol < 2. 



Ksea. = 0.156 




nconfigs 


= 198, nboot -- 


= 200, correlated 








(^minj^max) 


(9,15) 








(9,15) 








K1-K2 




\ 
/ 






3-V2(o|F'|l/) 




0.1585-0.1585 


0.0289 + 


11 
10 


0.0985 + 


25 
26 


1 fil 1 


82 
80 


0.4758 + 


99 
99 


0.1585-0.1580 


0.0314 + 


12 
10 


0.1008 + 


25 
27 


1 f\d7 ~^ 

VJ. ±\J^ i 


76 
78 


0.4680 + 


96 
99 


0.1580-0.1580 


0.0338 + 


12 
10 


0.1030 + 


27 
29 


1 (\R'^ ~^ 
VJ. xuoo 


67 
77 


0.4609 + 


93 
99 


0.1585-0.1575 


0.0337 + 


12 

10 


0.1029 + 


27 

28 


0.1685 + 


68 
77 


0.4603,^22' 


0.1580-0.1575 


0.0362 + 


13 
11 


0.1051 + 


28 
30 


0.1722 + 


73 
75 


0.4537 + 


94 
92 


0.1585-0.1570 


0.0361 + 


12 
11 


0.1049 + 


28 
28 


0.1725 + 


72 
74 


0.4526 + 


96 
95 


0.1575-0.1575 


0.0386 + 


13 
12 


0.1073 + 


28 
31 


0.1763 + 


69 
75 


0.4472 + 


90 
94 


0.1580-0.1570 


0.0385 + 


13 
12 


0.1072 + 


28 
30 


0.1763 + 


68 
74 


0.4465 + 


91 
95 


0.1575-0.1570 


0.0410 + 


14 
12 


0.1094;^ 


29 
■■il 


0.1806 + 


66 
71 


0.4405 + 


90 
92 


0.1585-0.1560 


0.0408 + 


13 
12 


0.1087 + 


27 
29 


0.1804 + 


71 
73 


0.4367 + 


96 
86 


0.1570-0.1570 


0.0434 + 


14 
14 


0.1116 + 


30 
32 


0.1849 + 


64 
69 


0.4343 + 


92 

88 


0.1580-0.1560 


0.0433 + 


14 
13 


0.1110 + 


29 
31 


0.1846 + 


63 
71 


0.4317 + 


92 
88 


0.1575-0.1560 


0.0458 + 


14 

13 


0.1133 + 


29 

32 


0.1891 + 


63 
66 


0.4266 + 


90 
89 


0.1570-0.1560 


0.0483 + 


15 
13 


0.1156 + 


29 
31 


0.1937 + 


62 
66 


0.4214 + 


89 
86 


0.1560-0.1560 


0.0534 + 


16 
15 


0.1199 + 


32 
31 


0.2030 + 


63 
68 


0.4103 + 


89 
90 
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TABLE X. Lattice results pseudoscalar and vector meson decay constants at Kgea = 0.1565. 



For all fits wo find 1 < x^/d.ol < 2. 



I^sea. = 0.1565 




nconfigs 


= 198, nboot = 


= 200, correlated 








(^minj^max) 


(9,15) 








(10,15) 






K1-K2 




\ 
/ 










p) 


0.1580-0.1580 


0.0317 + 


19 
17 


0.1024 + 


50 
43 


n 1 fi9 ^ 


12 
11 


0.481 + 


17 
18 


0.1580-0.1575 


0.0341 + 


20 
18 


0.1047 + 


48 
43 


yj. ±\jo 


11 
11 


0.472 + 


17 

18 


0.1575-0.1575 


0.0365 + 


21 
17 


0.1070 + 


50 
41 


n 1 fis ~^ 


11 
11 


0.463 + 


17 
18 


0.1580-0.1570 


0.0365 + 


21 

17 


0.1068 + 


50 

42 


0.169 + 


11 

11 


0.462 + 


16 
18 


0.1575-0.1570 


0.0390 + 


22 
18 


0.1091 + 


49 
40 


0.172 + 


12 
11 


0.454 + 


16 
17 


0.1580-0.1565 


0.0389 + 


22 
18 


0.1087 + 


50 
41 


0.172 + 


12 
11 


0.453 + 


16 
17 


0.1570-0.1570 


0.0414 + 


23 

19 


0.1113 + 


48 
41 


0.176 + 


11 
11 


0.446 + 


15 
17 


0.1575-0.1565 


0.0414 + 


23 
19 


0.1111 + 


48 
41 


0.176 + 


11 
11 


0.445 + 


15 
17 


0.1570-0.1565 


0.0439 + 


22 
19 


0.1133 + 


48 
:-!!) 


0.180 + 


9 
9 


0.438 + 


14 

16 


0.1580-0.1560 


0.0412 + 


22 
18 


0.1104 + 


50 
41 


0.176 + 


12 
11 


0.444 + 


15 
17 


0.1565-0.1565 


0.0463 + 


21 
20 


0.1153 + 


47 
39 


0.184 + 


9 
9 


0.431 + 


14 
15 


0.1575-0.1560 


0.0437 + 


23 
19 


0.1129 + 


49 
40 


0.180 + 


11 
10 


0.437 + 


14 
16 


0.1570-0.1560 


0.0462 + 


22 
19 


0.1151 + 


48 
38 


0.184 + 


9 
10 


0.430 + 


13 
15 


0.1565-0.1560 


0.0487 + 


22 
21 


0.1172 + 


47 
41 


0.188 + 


9 
10 


0.424 + 


12 
15 


0.1560-0.1560 


0.0512 + 


25 
21 


0.1191 + 


47 
41 


0.193 + 


9 
9 


0.417 + 


12 
14 
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TABLE XI. Lattice results pseudoscalar and vector meson decay constants at Kgea = 0.157. 

For all fits wo find 1 < x^/d.ol < 2. 



^sea — 0.157 




nconfigs 


= 198, nboot -- 


= 200, correlated 








(^minj^max) 


(10,15) 








(9,15) 








K1-K2 




\ 
/ 






3-V2(o|F'|l/) 




0.1575-0.1575 


0.0302 + 


10 
9 


0.0955 + 


27 
24 


0.1491 + 


65 
66 


0.4454+ 


79 

72 


0.1575-0.1570 


0.0326 + 


10 

9 


0.0979 + 


26 
24 


0.1538 + 


64 
70 


0.4388 + 


79 
75 


0.1570-0.1570 


0.0350 + 


11 
10 


0.1003 + 


26 
23 


0.1586 + 


62 
71 


0.4326 + 


75 
78 


0.1575-0.1565 


0.0349 + 


11 

10 


0.1001 + 


26 

23 


0.1584 + 


63 
72 


0.4317 + 


75 

78 


0.1570-0.1565 


0.0374 + 


12 
10 


0.1026 + 


25 
25 


0.1633 + 


67 
70 


0.4261 + 


72 
75 


0.1565-0.1565 


0.0398 1 


13 
11 


0.1050 + 


27 
27 


0.1680 + 


63 
66 


0.4201 + 


73 
73 


0.1575-0.1555 


0.0395 + 


12 
10 


0.1040 + 


27 
25 


0.1671 + 


67 
67 


0.4167 + 


72 
73 


0.1570-0.1555 


0.0420 + 


13 
12 


0.1067 + 


to to 

00 00 


0.1721 + 


66 
67 


0.4121 + 


72 
69 


0.15G5-0.1555 


0.0145 + 


14 
12 


0.1091;^ 


28 
28 


0.1770 + 


62 
(ii 


0.4071 + 


70 
(i7 


0.1555-0.1555 


0.0493 + 


14 
13 


0.1133 + 


26 
28 


0.1863 + 


63 
63 


0.3961 + 


68 
64 
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TABLE XII. Lattice results pseudoscalar and vector meson decay constants at 



K«,a = 0.1575. For all fits wo find 0.5 < x'^/(\.Q.i< 1. 



Ksea = 0.1575 




nconfigs 


= 198, nboot -- 


= 200, correlated 








(^minj^max) 


(9,15) 








(10,15) 












\ 
/ 






3-V2(o|F'|l/) 




0.1575-0.1575 


0.0251 + 


8 
9 


0.0894 + 


26 
28 


0.1286 + 


55 

61 


0.4443 + 


99 
99 


0.1575-0.1570 


0.0276 + 


8 
9 


0.0924 + 


26 
27 


0.1321 + 


57 
60 


0.4370 + 


99 

91 


0.1570-0.1570 


0.0301 + 


8 
9 


0.0953 + 


26 
26 


0.1361 + 


55 
55 


0.4305 + 


96 
86 


0.1575-0.1565 


0.0300 + 


8 
9 


0.0951 + 


27 

26 


0.1360 + 


53 

57 


0.4295 + 


95 
87 


0.1570-0.1565 


0.0325 + 


8 
9 


0.0979 + 


24 
26 


0.1403 + 


53 
57 


0.4239 + 


95 
82 


0.1565-0.1565 


0.0350 + 


8 
9 


0.1004 + 


23 
25 


0.1448 + 


53 
55 


0.4180 + 


91 
84 


0.1575-0.1555 


0.0347 + 


9 
10 


0.0994 + 


24 
27 


0.1439 + 


51 
57 


0.4141 + 


91 
84 


0.1570-0.1555 


0.0372 + 


8 
9 


0.1021 + 


23 
24 


0.1487 + 


52 
56 


0.4099 + 


90 
83 


0.15G5-0.1555 


0.0396 + 


8 
10 


o.ioig;^ 


23 
22 


0.1535 + 


52 

r)8 


0.4052 + 


88 
«5 


0.1555-0.1555 


0.0444 + 


10 
10 


0.1089 + 


21 
22 


0.1629 + 


47 
57 


0.3948 + 


83 
80 



TABLE XIII. Fit results for pseudoscalar, vector, nucleon and delta particles (in lattice 
units) from "symmetric" fits. 



particle 


a 


h 




xVd-O.f 


PS 


-12.407(380) 


1.9666(540) 




1.2/2 


particle 




c 


e 


xVd.o.f 


V 


0.3300(90) 


4.070(250) 





1.5/2 


V 


0.2928(410) 


6.44(2.40) 


-33.28(34.1) 


0.27/1 


N 


0.5012(190) 


6.960(460) 





2.7/2 


N 


0.4246(750) 


11.84(4.8) 


-67.99(69.1) 


1.5/1 


A 


0.5851(240) 


6.482(670) 





6/2 


A 


0.444(110) 


15.16(6.3) 


-118.9(82.0) 


3.5/1 




0.0496(34) 


0.888(101) 





1.3/2 


u 


0.0423(137) 


1.354(778) 


-6.58(11.15) 


0.89/1 


v/p 


0.302(110) 


-0.372(290) 





0.37/2 


v/p 


0.294(43) 


0.170(2.204) 


-7.45(31.44) 


0.31/1 
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TABLE XIV. Values of the coupling constant. We do not take into account the light-quark 
dependence of the plaquette in the determination of the strong coupling constant. 







0^(3.41/0) a-jjgiTr/a) 




P = 5.6, Nf = 2 


(0.43012,0.42927,0.42837,0.42749) 


0.167 0.150 


0.215 


/5 = 6.0, A^i. = 


0.406318 


0.152 0.138 


0.205 


TABLE XV. 
sector, according 


Fit results for the masses of pseudoscalar and vector particles in 
to eqs. 15, 17 and 18 


the strange 


PS (linear fit) 


o! 


h' 


xVd.o.f 




3.93(12) 1,01(11) 


26/31 


V(linear fit) 


^/crit ^ 


d' 


xVd.o.f 




0.3300(93) 4.07(25) 0.948(31) 


13/31 


V(quadr. fit) 




d' 


xVd.o.f 




0.2928(412) 6.44(2.4) 0.908(145) 


3/29 




/' 


9' 




-33.28(34.3) -2.19(4.79) 2.25(2.03) 


TABLE XVL Collection of results for «;Strange_ 


particle 


strange 




strange 
^quad 


K 


0.15654(11) 




0.15694(42) 


K* 


0.15561(14) 




0.15590(57) 


4> 


0.15563(14) 




0.15598(50) 
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TABLE XVII. Lattice results in the strange quark sector. 



UUocl VcLUic 


iin TIT T r~i 'Kvi 
llli liL LU Itip 




mJ-dU- ill LU iitp 








n 9zL1 










TYi Tj'^ I /Ti 1 

I'l'K* W) 


U.OOOl ± -L ( 












OA^fOR^ 








U.Oc/Ot fJU ) 


mK*{K*) 


0.288(11) 




0.346(48) 


m,K{K) 


0.215(6) 




0.192(26) 


m^{K) 


0.407(13) 




0.356(53) 


mK'iK) 


0.371(11) 




0.326(47) 


fK{(t>) 


0.0633(33) 




0.0615(38) 


-I 1 r / 1 \ 


0.2742(93) 




0.2777(97) 




0.0633(33) 




0.0619(38) 




0.2741(93) 




0.2769(104) 




0.0591(33) 




0.0572(34) 




0.2832(97) 




0.2871(105) 


TABLE XVIII. 


Values of the inverse lattice spacing 


obtained from different observables 


(at the hght quark 


mass). The lattice value of mjv results from a quadratic extrapolation, the 


value of /tt stems from a linear extrapolation. 






observable 


m(4ie*) 




a-i 


mp (linear fit) 


0.334(9) 




2.30(6) 


Viip (quadratic fit) 


0.297(41) 




2.58(37) 


mjv (rho linear) 


0.435(72) 




2.16(40) 


(rlio liiioar) 


0.0505(31) 




2.62(18) 


TABLE XIX. Physical results in the light quark sector. 


mN[GeV] 


mA[GeV] 


fAGeV] 




linear vector 


1.00(16) 


1.05(23) 


0.116(8) 


0.302(11) 


quadr. vector 


1.12(25) 


1.17(28) 


0.130(21) 


0.302(11) 
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TABLE XX. Physical results for masses in the strange quark sector. 



mK{K*)[GeV] 


mKmOeV] 


mK*{K)[GeV] 


mK*{(l))[GeV] 


m4,{K)[GeV] 


m^{K*)[GeV] 


linear vector 


0.596(9) 


0.595(9) 


0.853(3) 


0.894 


0.937(6) 


1.020 


quadr. vector 


0.633(47) 


0.623(44) 


0.842(12) 


0.890(4) 


0.918(24) 


1.029(9) 



TABLE XXI. Physical results for decay constants in the strange quark sector. 



fKip)\GeV] 






l/.fJ.o) 




l//o(A') 


linear vector 


0.1456(71) 


0.1457(71) 


0.1360(72) 


0.2742(93) 


0.2741(93) 


0.2832(97) 



TABLE XXII. Physical results in the quenched sector. Fits to nucleon and A are quadratic, 
all other fits are linear. To set used the Kaon mass for fx and uik*, and the 

mass for !//</,, rriK and J. 



m,N[GeV] 


m,A[GeV] 


fn[GeV] 






1.0(31(67) 


1.301 ((33) 


0.132.5(18) 


0.3271 (.52) 




niKlGeV] 


niK* [GeV] 


MGeV] 




J 


0.5590(64) 


0.8647(23) 


0.1518(46) 


0.2982(47) 


0.38(1) 
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FIGURES 



Pseudoscalar(ls) 



0.2 0.4 



0.2 0.4 



Pseudoscalar(ss) 



0.2 0.4 



Vector(ls) 



0.2 0.4 



Vector(ss) 




FIG. 1. The standard error of the masses of pseudoscalar, vector and nucleon as a function 
of the inverse block size 
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0.9 - 
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0.7 - 



0.6 - 



0.5 - 



0.4 - 



0.3 
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N 







■'■ * -■■ ^ ? >^'- J: 



V 



X 



>^ 4 I 



PS 



^ ^ ^ 



8 10 12 14 



t 



FIG. 2. Effective masses at Kg^a = 0.1575. * are smeared-smeared and x are smeared-local 
data. The results of our fits to the smeared-smeared correlators (not to the effective masses!) 
are indicated by solid lines. 
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0.01 0.02 0.03 0.04 0.05 0.06 

"^sea 

FIG. 4. rriN and 

771a as a function of Ulsea 

(in lattice units). 
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1^sea=0- 560 ^ 

K3ea=0.1565 - ■ 

Ksea=0;570 - 

Ksea=0-1575 - ♦ 




0.05 0.1 0.15 0.2 

FIG. 7. Mass of the vector meson as a function of mi. 
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1.6 - 



1.4 



1.2 - 



0.8 



0.6 - 




FIG. 8. Comparison of Np = 2 (circles) and quenched results (squares) in the light sector 
with experiment. The data is normalized to its experimental values, namely Mn = 938MeV, 
Ma = 1232MeV, U = 132MeV, l//p = 0.199a/2. To set the scale we used the linear fit to the 
vector meson trajectory. 
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1.6 - 



1.4 



1.2 - 



0.8 



0.6 - 




FIG. 9. Comparison oi Np = 2 (circles) and quenched results (squares) in the strange sector 
with experiment. The data is normalized to its experimental values, namely Mk = 495MeV, 
Mk* = 892McV, fx = leOMeV, 1//$ = 0.234 and J = 0.48. To set the scale we used the 
linear fit to the vector meson trajectory. 
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